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One of the most spectacular experimental findings in the fractional quantum Hall effect is evidence for an 
emergent Fermi surface when the electron density is nearly half the density of magnetic flux quanta {u = 1/2). 
The seminal work of Halperin, Lee, and Read (HLR) first predicted that aiv = 1/2 composite fermions—hound 
states of an electron and a pair of vortices—experience zero net magnetic field and can form a “composite Fermi 
liquid” with an emergent Fermi surface. In this paper we use infinite cylinder DMRG to provide compelling 
numerical evidence for the existence of a Fermi sea of composite fermions for realistic interactions between 
electrons at = 1/2. Moreover, we show that the state is particle-hole symmetric, in contrast to the construction 
of HLR. Instead, our findings are consistent if the composite fermions are massless Dirac particles, at finite 
density, similar to the surface state of a 3D topological insulator. Exploiting this analogy we devise a numerical 
test and successfully observe the suppression of 2k f backscattering characteristic of Dirac particles. 


Electrons confined to a two-dimensional (2D) plane in a 
strong magnetic field organize into a plethora of remarkable 
phases in which electron correlations play the key role. A 
global understanding of these phases is provided by topologi¬ 
cal excitations called composite fermions (CF).'’^ One of the 
most intriguing phases occurs in the half-filled {v — Lan¬ 
dau level, in which composite fermions were predicted by 
Halperin, Lee, and Read (HLR) to form a metallic state. ^ 
Experiments corroborated this understanding by identifying 
signs of a Fermi surface in spite of the intense magnetic 
field."^^ Despite the tremendous success of the HLR theory, 
one aspect of the v = \ phase has remained an enigma: 
a particle-hole (PH) symmetry of the quantum Hall prob¬ 
lem projected into the spin-polarized lowest Landau level 
(LLL).*^'^ The HLR theory attaches fluxes to electrons rather 
than holes, and is not confined to the LLL, breaking PH sym¬ 
metry. One possibility is that PH symmetry is spontaneously 
broken,'"^ as is believed to occur in the excited Landau level 
(e.g. GaAs at filling v = |),*^“** leading to the Moore-Read 
phase.'® 

After two decades of study, theorists have proposed a 
radical twist to this picture, forging deep connections with 
the physics of three-dimensional (3D) topological insulators 
In the new picture due to Son,^" composite fermions 
are massless Dirac particles, similar to electrons on the surface 
of a TI. The CFs are coupled to an emergent gauge field, so 
the new theory retains the non-Fermi-liquid aspects of HLR— 
in fact, this “Dirac-CFL” is equivalent to finite-density QED 3 . 
In a TI, the masslessness of the Dirac fermions is protected 
by the Kramers time-reversal symmetry. Here time-reversal 
is absent, but the PH symmetry that trades occupied and un¬ 
occupied electronic states at jz = ^ plays an equivalent role. 
It is proposed that when PH acts on the CFs, it behaves ex¬ 
actly like Kramers time-reversal symmetry, exchanging states 
on opposite points of the Fermi surface and forbidding a CL 


mass term.^" Rather than an enigma, PH symmetry now plays 
a starring role. 

We address the following two questions. First, is the half- 
filled Landau level consistent with a composite Fermi liquid, 
either in its original form or the Dirac revision? Second, is 
PH symmetry preserved, and what are the measurable effects 
if the composite fermions are indeed Dirac particles? Using 
large scale density matrix renormalization group (DMRG) nu¬ 
merical simulations,^^ we provide strong evidence for the for¬ 
mation of a CFL in the LLL with realistic Coulomb inter¬ 
actions. We observe a single, nondegenerate Fermi surface 
consistent with a Luttinger count of CFs, and signatures of an 
emergent gauge field. In the idealized limit of a single Landau 
level, particle-hole symmetry is preserved. Combining these 
numerical observations with new theoretical insights into the 
microscopic action of PH, we argue this is logically suffi¬ 
cient to conclude that the CL is a Dirac fermion. However, to 
demonstrate this point directly, we test for the analog of a clas¬ 
sic signature of the TI surface—suppression of 2kF backscat¬ 
tering off impurities that preserve time-reversal symmetry. 

We indeed find that backscattering off PH-symmetric poten¬ 
tials is suppressed in the CFL, reappearing only when PH- 
breaking perturbations are introduced. We also find that the 
effect of the small PH-breaking perturbations present in ex¬ 
periments is very weak, suggesting that PH symmetry is ex¬ 
perimentally relevant at z/ = ^. 


I. MODEL AND METHODS 

We study electrons on an infinitely long cylinder, in an 
external magnetic field B, interacting via Coulomb repul¬ 
sion. There are several advantages to the infinite cylinder ge¬ 
ometry: there are no edge effects, PH symmetry can break 
spontaneously, and there can be algebraic correlations. For 
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our numerical simulations,^® we truncate the interaction as 
^ and project into the lowest Landau level. 
We assume the system is spin-polarized, which is observed 
to occur in the high field regime.We fix A = 6£b, which 
is large enough to capture the Coulomb interaction physics 
but small enough to avoid wrapping effects on a cylinder. 
Throughout, the magnetic length £b = \/hf{eB) is set to 
1 ; thus lengths are given in multiples of momenta are 

measured in units of We use the coordinate x along the 
infinite length and y around the finite circumference Ly of the 
cylinder. 

The states we study are critical and have algebraic corre¬ 
lations along the cylinder, as well as infinite bipartite entan¬ 
glement. However, infinite-DMRG always introduces a finite 
length scale to the system due to the finite bond-dimension 
X used in the DMRG variational ansatz. Instead of the usual 
“finite size scaling,” we perform the so-called “finite entangle¬ 
ment scaling” by running DMRG at various bond dimensions 
X « 10^-10'^, and extrapolating the results as x 
This way we can analyze the ^ ^ oo limit and extract the 
critical properties of the system on the infinite cylinder. 

We describe further details of our DMRG setup in Appen¬ 
dices B and C. Here, we only note that the momentum Ky 
around the cylinder is conserved in the DMRG, so it must be 
chosen correctly in order to find the ground state. We find 
that the momentum Ky of the ground state depends on the 
circumference Ly. Only two such momenta appear, and for 
reasons explained in App. C we label these momentum sec¬ 
tors by “0110” and “0101.” Interestingly, as we show below, 
these sectors allow us to access different boundary conditions 
for the CFs, even though the CFs are emergent fields. 


II. MAPPING THE FERMI SURFACE OF EMERGENT 
FERMIONS 

In the 2D limit, the composite fermions are expected to 
form a circular Fermi surface with kp — ^ 3.1 |, which 

we can study with our numerics. While the cylinder is infinite 
in X, the finite Ly quantizes the momenta ky into a discrete 
grid with spacing jLy. Thus, instead of a circular Fermi 
surface, we expect a set of “wires” labeled by ky, illustrated 
in Fig. 2(a). As Ly is increased, the available ky’s get closer 
together and the number of wires increases, better approx¬ 
imating the 2D system. 

We get a sense of how well we approach the 2D limit in 
Fig. 1, where we plot the electron density-density correla¬ 
tions D{q} = ( :6pqSp-q :) as a function of momenta q at 
Ly = 24. [See App. C 2 for details on measuring D(q).] The 
density structure factor in the 2D limit is expected to be ro- 
tationally symmetric with a singularity on a circle of radius 
2kF. While we only have access to data at discrete qy, our 
measurements still approximate a circular shape, with distinct 
singular features near |q| = 2kF. One expects a singularity 
^(b) |lHl~2fcF| with a = I for free fermions and for the 
CFL with Coulomb interaction, but a is modified for the CFL 
with short-range interactions.Future studies approaching 
closer to the 2D regime should allow such predictions to be 
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FIG. 1. Density-density correlations D(q) = (: SpqSp-q :) plotted 
as a function of 2D momenta q = {q^;, qy). The cylinder circum¬ 
ference is Ly — 24, where we see eight slices through the Fermi 
sea. Descendants of a singular 2kF circle in D(q) show up in cuts 
at qy between —7(^) and 7(^) (see also Fig. 7). At this large 
circumference the correlations approach those of a 2D system. 


tested quantitatively. 

We determine the number and lengths of the CF wires 
from £l(q). The electron density will generically couple to 
CF scattering processes (allowed by the symmetries), i.e., 
^k^CF;kV^CF;k-i-q + • ■ •, where ellipses denote 
higher-body processes. A transition of momentum q across 
the Fermi surface will contribute a singularity to the structure 
factor D{q). Since k and k -F q are restricted to the wires, 
the singularities in D{q) can then be used to determine the 
configuration of the wires. For example, the singularities in 
D{qx: Qy = 0) contain transitions within the wires (i.e., fixed 
ky), and reveal the lengths of the wires inside the Fermi sea. 
At other qy, the singularities correspond to processes which 
connect the ends of the wires whose ky momenta differ by qy. 

Slices of the density structure factor for Ly = 13 are 
shown in the lower panel of Fig. 2(b). Note that here we 
plot D{q) = e'l /^I7(q), which has the same singularities 
as I7(q) but is more conveniently scaled (see App. C2). Vi¬ 
sual inspection reveals some of the expected singularities, and 
we can increase the contrast by taking a “fractional” deriva¬ 
tive with respect to q^.^^ We calculate an p* order derivative 
by multiplying the real-space correlations by \x — x'\'^ before 
Fourier transforming. In the upper panel we show the results 
for various p G (0.5,1.5), with rj chosen for each singularity 
individually. This method reveals many singularities which 
are barely visible in the raw data. 

Interestingly, while the physical electrons have periodic 
boundary conditions, the composite fermions need not. An 
arbitrary flux $int of the emergent internal gauge field can 
thread the cylinder, quantizing the CF momenta as ky G 
Being a dynamical degree of freedom, the 
emergent flux $int will adjust to lower the energy of the filled 
Fermi sea, depending on Ly. Only two cases respect 180° ro¬ 
tational symmetry: periodic boundary condition (PBC) with 
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FIG. 2. (a) Boundary conditions on the composite fermions. Our system can be described by a number of “wires”—slices through the 2D Fermi 
sea—at various fixed ky. The number of wires is dictated by Ly and the boundary condition (BC) of the composite fermions. An odd number 
of wires (upper panel) corresponds to PBC (periodic); an even number of wires (lower panel) corresponds to APBC (antiperiodic). (b) Mapping 
the Fermi surface via the structure factor. The lower panel shows the LLL density-density correlations e'’ ^^D(q) = { :(5pq5p_q:) 

measured on a cylinder with Ly = 13. The singularities arise from CF scattering processes across the Fermi surface. The observed scatterings 
are illustrated in the inset, with colors corresponding to qy. The upper panel shows the derivatives of the correlator, which aids in determining 
the location of the singularities, (c) Testing Luttinger’s Theorem. Luttinger’s theorem states that the area enclosed by the Fermi surface 
is related to the particle density. On a cylinder the “area” is given by the sum of the length of each wire in momentum space, which we 
determined from singularities in plots like Fig. 2(b). We define Qm to be the length of the Fermi sea slice at ky = ^m, and plot the resulting 
sums for various circumferences against the Luttinger’s prediction. Note that we use the relation Q-m = Qm, a consequence of rotation 
symmetry. There is excellent agreement between our data and the Luttinger count. 


int = 0 and anti-periodic boundary condition (APBC) with 
As shown in Fig. 2(a), the boundary condition dic¬ 
tates the parity of the number of wires: PBC yields an odd 
Nyj, while APBC yields an even Nyj. 

The observed singularities at Ly = 13 are in exceptional 
agreement with an Nyj = 4 model with APBC. The inset of 
Fig. 2(b) shows arrows for the expected singularities, all of 
which are observed in the main plot. We obtain similar data 
for other circumferences Ly (see App. D), with Fermi surfaces 
cut by Nyj = 4, 5, 6 , and 8 wires.We find the parity of 
Nyj depends only on the momentum sector used to initialize 
DMRG: 0101 always yields odd, while 0110 always gives 
an even N^j. This correspondence between the momentum 
sector and the BC is consistent with the CF theory, where it 
arises from an anomaly in finite-density QED 3 (see App. G). 

A useful check on our conclusions so far is a comparison to 
Luttinger’s theorem^® for the composite fermions. The elec¬ 
tron density per unit-length of the cylinder is piu — Ly .^^^^ . 
The Luttinger count requires that 27rpiD be equal to the vol¬ 
ume of the Lermi sea (here the sum of the lengths of the 
wires), as predicted by both the HLR and Dirac-CLL theories 
at 1 / = i. We find excellent agreement with this prediction in 
Lig. 2(c). Our findings are in contrast with a recent sugges¬ 
tion in Ref. 40 that the Luttinger count is violated by the CL 
model-wavefunctions. 

We can learn more about the composite fermions by mea¬ 
suring the central charge c. This is determined from the 
scaling between the entanglement entropy S and the DMRG 
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FIG. 3. Entanglement entropy vs. the correlation length at a variety 
of circumferences Ly. Different data points at the same size corre¬ 
spond to different bond dimensions x = 600-12000. Both S and ^ 
increase as x is increased. For a quasi-ID critical system, we expect 
a linear relationship between S and log ^ with the slope proportional 
to the central charge c [see Eq. (1)]. The dashed lines (from bottom 
to top) in the figure correspond to c = 3, 4, 5, and 7. 

ansatz correlation length f (the effective cutoff length) as the 
bond dimension x is increased. While in a gapped state S and 
f converge with x, for ^ critical state both quantities diverge. 
The central charge is extracted using the relation"^' 

S'=-log ^-f const. (1) 
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Figure 3 shows S vs. log ^ as x is varied. The dashed lines cor¬ 
respond to c of 3, 4, 5, and 7. Clearly c increases with Ly, with 
each new wire adding 1 to the central charge. The preceding 
analysis found = 4, 5, 6 , and 8 for these systems respec¬ 
tively, which allows us to deduce the relation c = N^j — 1. Our 
data rules out an ordinary Fermi liquid, for which c = iV^,. In¬ 
stead, it confirms the emergence of a gauge field in the com¬ 
posite fermion liquid; in the quasi-ID limit, the effect of the 
gauge field is to gap out the total gauge charge mode, reducing 
c by one (see App. I). 


III. PARTICLE-HOLE SYMMETRY AND THE ABSENCE 
OE 2kF BACKSCATTERING 

Our numerical findings thus far are expected of both the 
HLR and Dirac-CFL phases, so to distinguish them we turn to 
particle-hole symmetry. The symmetry is generated by TH, an 
anti-unitary operator which swaps creation operators with 
annihilation operators c: 

| 000 ---) ^ ( 2 ) 

Particle-hole is an exact symmetry of the Hamiltonian when 
the interaction is projected into the LLL, which is justified 
when the cyclotron energy hujc —oo. In reality PH is weakly 
broken by the finite ratio of the Coulomb energy to the cy¬ 
clotron energy. In App. E, we provide quantitative evidence 
that this breaking is weak in typical GaAs samples. 

The HLR constraction breaks PH, while the Dirac-CFL ex¬ 
plicitly preserves PH. Our first test is for spontaneous PH sym¬ 
metry breaking. We note that PH can be spontaneously broken 
in our infinite cylinder numerics, as it is a discrete symme¬ 
try that we do not explicitly enforce. Indeed, our numerical 
DMRG runs clearly break PH at v — |, randomly selecting 
either the Pfaffian or anti-Pfaffian state.Lurthermore, our 
numerics will never produce a symmetric superposition (“cat 
state”) of the HLR and its PH-conjugate (see App. C). 

To test for spontaneous symmetry breaking we compute 
the overlap between the ground state and its PH conjugate, 
('I'|7^ I'F) = (1 — which should decrease exponen¬ 

tially with the number of orbitals (i.e., system size) Ajjrb- 
DMRG directly computes the overlap per orbital 1 — e; non¬ 
zero e > 0 indicates symmetry breaking. Lor comparison, we 
calculate the same overlap in the n = 1 Landau level (i.e., 
iz = I mentioned above).As an example, at Ly = 16 for 
z/ = I we find e < 6 X 10 while for v = | we find 
e « 0.022(2). At 1 / = I this implies that the total overlap on a 
16 X 16 toms with 7Vo,b — 16^/(27r) is 0.978'^"'’ « 0.4, so the 
PH-breaking is in fact quite strong. The difference between 
V = \ and z/ = I is also manifest in qualitatively different 
“entanglement spectra,”"^^ as illustrated in Pig. 4. Similar re¬ 
sults hold at Ly = 13,18. Therefore, barring some transition 
at even larger Ly, we conclude that the v = \ CPL state is 
PH-symmetric. 

Given the presence of PH symmetry, we can ask how the 
CPs transform under this symmetry. According to Refs. 20- 
22, PH acts like time-reversal symmetry on the CPs implying 



FIG. 4. Contrasting orbital entanglement spectra of the n = 0 and 
n = 1 Landau levels at half-filling (corresponding to experiments 
at z/ = 1 and z^ = | respectively). Data is at Ly = 19; n = 0 
is the gapless CFL phase, and n = 1 is consistent with the gapped 
Pfaffian phase. Lent corresponds to the eigenvalues of the reduced 
density matrix for the left half of the cylinder. Each eigenvalue is 
associated with a momentum around the cylinder Ky, just like the 
energy spectrum of a disk. PH acts as a reflection Ky -H- —Ky. The 
characteristic chiral “dispersion” at n = 1 clearly breaks PH, while 
n = 0 is PH-symmetric. 

a twofold Kramers degeneracy whenever an odd number of 
CLs is present. Indeed, this feature is already apparent by 
the very definition of VH, since a system with A(,jb orbitals 
transforms as (cf App. A) 

Thus at filling i/ = i, adding a pair of fluxes plus an 
electron—i.e., a composite fermion—changes by —1. 

The Kramers degeneracy of the CLs leads to their twofold 
“pseudospin” degree of freedom, which are exchanged via 
PH. Two kinds of low energy gapless theories are possible. 
The first is a pair of CF Fermi surfaces which are exchanged 
by PH (or possibly Rashba split). The other is a single nonde¬ 
generate Fermi surface in which the CF pseudospin is locked 
to the momentum, that is, a Dirac cone. Our numerical find¬ 
ing of a single Fermi surface with the correct Luttinger vol¬ 
ume implies the latter proposal is realized, consistent with 
Refs. 20-22. 

One of the most striking consequences of a Dirac cone is 
that composite fermions accumulate a tt Berry phase when 
circling the Dirac point. In topological insulators, this Berry 
phase prevents precise 2 fci? backscattering off time-reversal- 
symmetric impurities on the surface. If our system is a Dirac 
CFL, we should observe the same effect, with T’H playing 
the role of time-reversal operator. At the level of equal-time 
correlation functions, we expect that a Hermitian, PH-even 
operator P will not have strong singularities in (PqP_q) at 
|q| = 2kp in 2D (cf. App. I). The electron density studied 
earlier is odd under PH, so does not have this property. A can¬ 
didate PH-even observable is P(r) = Sp{r). A contri¬ 
bution to this correlation function (see App. D for details) is 
shown in Fig. 5, measured for the same four-wire APBC CFL 
as in Fig. 2(b). The measurements are done at qy — 2tt/ Ly, 
so as to probe the CF scattering from the right Fermi point at 
ky — tt/Lj, to the left point at ky = —ir/Ly, which corre¬ 
sponds to exact 2k F backscattering in 2D. 

Despite showing many other singularities, all of which can 
be accounted for by various multiple-CF scattering processes 
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FIG. 5. The extinction of 2k f backscattering off PH-symmetric im¬ 
purities. A QH bilayer with chemical potential imbalance fj, allows 
us to continuously tune from a PH-symmetric model (/r = oo) to a 
PH-broken one {p finite). We compute the correlation function of 
a PH-even operator {PqP_q) for qy = ^ al Ly = 13 and plot 
its derivative with respect to to bring out singularities. At the 
PH-symmetric point, there are many singularities in q^ (see App. D 
for analysis), but noticeably absent is any kink at |q| = 2kF- This 
demonstrates the Dirac structure of the CF Fermi surface: P is even 
under PH, while scattering a CF across the Fermi surface to its an¬ 
tipode is PH-odd in the Dirac theory. At finite chemical potential ^ 
the bilayer setup explicitly breaks PH symmetry, and a kink at 2k f 
continuously reappears. 


(see App. D), the result is perfectly smooth in the vicinity of 
the (already measured) momentum corresponding to this exact 
2kF backscattering. We also confirm the absence of the exact 
backscattering in the PH-symmetric model for a wider cylin¬ 
der with Ly = 16 which realizes five-wire CFL (see App. D). 

If the absence of the backscattering is truly due to the 
PH symmetry, rather than some peculiarity of P, we expect 
that the backscattering should return if PH symmetry is ex¬ 
plicitly broken. We break PH symmetry by adding a second 
quantum well at a distance 1£g above the first, with Coulomb 
intra- and inter-well interactions. The second well has a chem¬ 
ical potential p relative to the first one, and electrons can tun¬ 
nel between the wells. When p = oo, electrons remain in the 
first well, and we recover PH symmetry. As p decreases, elec¬ 
trons tunnel to the second well, breaking PH symmetry, which 
should induce backscattering. 

In Fig. 5 we show data confirming this hypothesis for Ly = 
13. The tunneling strength is fixed at t = 0.01 in units of 

/{AireiB)- The system remains in an effectively one-well 
CFL phase for p > 0.05 [for smaller p the system becomes a 
Halperin (331) state ]. As pis decreased within the CFL phase, 
the 2kp singularity reappears. Note that the measured central 
charge of the CFL phase remains unchanged; thus, the gapless 
Fermi surface is stable against PH-breaking perturbations (see 
also App. I). 

Our findings appear to violate the famous 2D fermion dou¬ 
bling theorem: a single time-reversal-symmetric Dirac cone 
is anomalous, so cannot be realized in 2D. However, as noted 


in Ref. 22 our composite Dirac fermion is coupled to an 
emergent gauge-field with unusual compactification condi¬ 
tion, which “cures” the anomaly. The half-filled LLL has 
yet another anomalous property. Just as the fully symmet¬ 
ric surface of a topological insulator must be nontrivial—that 
is either gapless or topologically ordered—a PH-symmetric 

state of the half-filled Landau level must be nontrivial. In- 

2 

deed, it must have a Hall conductance of and if the 

state is gapped, this requires fractionally charged excitations 
and hence topological order. In fact PH here behaves exactly 
like time-reversal symmetry on the surface of a 3D topologi¬ 
cal phase (in class AIII—see App. H). How does this occur in 
a purely 2D setting? For any symmetry that is locally imple¬ 
mented, one can always obtain a symmetric and trivial product 
state. The key observation resolving this apparent paradox is 
that PH symmetry of the LLL is special, in that its action is 
nonlocal. The nonlocality is ultimately tied to the fact that 
the Landau level orbitals ipj (r) cannot be localized due to the 
topological nature of the LLL. 

We conclude with a number of open questions of relevance 
to both experiment and theory. First, given the experimen¬ 
tal success in observing the phenomenology of a Dirac cone 
in TI surfaces and a CF Fermi surface in GaAs, it would be 
extremely interesting to find experimental probes of PH sym¬ 
metry and the potential Dirac nature of CFs. Existing experi¬ 
ments are already of some relevance, such as recent measure¬ 
ments in Ref. 43 of CF “geometric resonances” induced by 
small deviations of H-field away from half-filling. DMRG 
could be of use in guiding and interpreting such experiments, 
for instance by computing static structure factors, impurity re¬ 
sponses, and the behavior of the CFL atv = ^ + 5. 

Second, Son has proposed a PH-symmetric version of a 
paired phase, the “PH-Pfaffian,”^° which has previously been 
proposed as the surface topological order of a PH-symmetric 
(class AIII) 3D topological superconductor.^"^ While our re¬ 
sults appear to rule out this possibility in the n = 1 LL of 
GaAs, it would be interesting to search for such a phase in 
broader phase diagrams of PH-symmetric models. 

Last but not least, similar theories with a surface of emer¬ 
gent gapless fermions coupled to an emergent gauge field arise 
for other exotic phases with itinerant fractionalized excita¬ 
tions, such as spin liquids with a spinon Fermi sea, Bose- 
metals, and electron non-Fermi-liquid metals. Much recent 
theoretical effort has aimed to clarify the status of such field 
theories,though it remains not fully settled away from 
artificially controlled limits. Unbiased numerical studies of 
the CFL thus bear directly on open questions for all these 
other non-Fermi liquids. Recent numerical studies"^^ explored 
quasi-lD ladder descendants of various non-Fermi liquids. 
Thanks to many innovations in the DMRG for FQHE, the 
present CEL work goes to effectively much wider strips and 
is much closer to the 2D physics than the previous studies. It 
would be useful to push the numerical CEL study yet closer to 
2D and develop scaling analysis tools for addressing 2D ques¬ 
tions, such as detailed characterization of the 2kF singularity 
in the structure factor. Time-dependent DMRG could poten¬ 
tially study the dynamical properties of a non-Fermi liquid, 
which has not yet been investigated numerically. 
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VH? = —1 relation for CFs (unpublished). 
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Appendix A: Particle-hole symmetry and Kramers degeneracy of composite fermions 

Here we discuss important formal properties of the particle-hole symmetry. One can define the following particle-hole (PH) 
transformation; 

c(r) —>• c’l’(r) , 

TH : c1'(r) ^ c(r) , (Al) 

j —>■ —i, 

where c(r) is the electron destruction operator in the continuum. Note that VH is an anti-unitary transformation. Now consider 
restricting to the lowest Landau level. Then, if cj creates a particle in the lowest Landau level orbital with wavefunction 
thus 

c] = ^ (r) (r) , ^ 
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we have that 


TH : Cj Cj , i ^ —i . (A3) 

Note that even though c/s depend on the choice of orbitals, the PH transformation is basis-independent. 

We also need to specify the action of particle-hole transformation on the vacuum state |0), which is defined as cj |0) = 0 for 
all j (i.e., the empty state). Let us denote |0') = TH |0). Under particle-hole the appropriate condition is cj |0') = 0 for all j, 

which implies that the transformed state is completely filled, |0') = Cj |0). In the latter equation, we understand some fixed 
ordering of fermion fields, e.g., 

|0')=7^|0) = c1c|...4J0) , (A4) 


where Norh is the total number of orbitals. 

The key result of this section, which we will show momentarily, is 


_ ^_J^^Aforb(JVorb-l)/2 


(A5) 


As a warm-up example, consider a simple case of one electron with two orbitals. Call the two orbitals ci, C 2 , and label the 
states |1) = c| |0) and |2) = c\ |0). Now, under particle-hole we have 

|l) = cl|0) 7^|1)=ci(cI4|0)) = + |2), 

|2) = 4|0) ^ m\ 2 )=c 2 iclcl\ 0 )) = -\l) . 


Thus TH^ = —1 for A),rb = 2. Since VH is an anti-unitary operation, there is no way to gauge out the minus sign, and this 
implies a Kramers degeneracy of the two states. 

This can be readily extended to an arbitrary number of orbitals. Consider a general basis state that can be represented as: 


•^orb 

I'l')=n io). (A6) 

i=i 

where nj G {0,1}. Acting twice with particle-hole yields for some number p. We first argue that p is 

independent of the number of creation operators Since T’H^CjT’H~‘^ = c|, one can commute the TH^ operator 

across all the creation operators: 


^orb ^orb 

|4') =m^Y[ lo) = n lo) = 1^) (o|^» ■ (A7) 

i=i i=i 

Therefore |0) picks up the same sign 77 under as that for every state with orbitals. To compute 77 , see that 77 10) = 
TH^ |0) = C 1 C 2 • • • CAr„,^c|c 2 • • • |0). Exchanging the raising/lowering operators until cj is next to takes j — 1 swap 

operations, for a total of ~ 1) = A),rb(A^orb ~ l)/2 negative signs. Therefore, we have shown that 77 = (—l)^orb(tVorb-i )/2 

as advertised in Eq. (A5). Again, we emphasize that the sign of is determined by A{,jb and not by the number of electrons. 
The sign is reversed each time we add two orbitals (i.e., two flux quanta). 

Particle-hole takes a state at filling ly, where 0 < 17 < 1, and transforms it into a state at filling \ — v. A state can only be 
PH-symmetric if it is at half filling {u = 1/2); thus Worb = 2Ne where Ne is the total number of electrons. In such case the sign 
depends only on whether Nf. is even or odd: (_i) 2 Afe( 2 Afe-i )/2 _ Pqj. particle-hole-symmetric states, we claim that 

states with an odd number of electrons must be doubly degenerate, which can then be interpreted as arising from a Kramers-like 
degeneracy of composite fermions. This may appear strange since time-reversal symmetry is broken by the magnetic field, and 
we always consider just one spin component of the electrons. However, one can define T' = VH as a new kind of time-reversal, 
which acts like regular time-reversal symmetry on the composite fermions leading to the Kramers degeneracy. 

Given Eq. (A5), one may ponder on the cases when A/ib is odd. If we naively start with two orbitals with VH^ = — 1 and 
take the square root, we may be mislead into thinking that VH^ = ±7 for a single orbital. Of course, this is false: Eq. (A5) 
gives VH^ = 1 in this case, which we can also verify directly. Indeed, VH |0) = e*“ |1) = e*“c^ |0) implies that VH^ |0) = 
VH{e^°‘V |0)) = e~^°'c{VH |0)) = |0) and thus VH^ = 1 on a single orbital. The resolution to this seeming paradox is that 
when acting on a single orbital TH is a fermionic operator—changing the fermion parity number—and thus does not have a 
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trivial product structure when combining orbitals. Labeling the two orbitals A, B, and defining operator TH = T’Ha'PHb, we 
obtain 




{THaTHb){THaPHb) = -VHaVHaVHbVHb = -TH\VHI = -1 , 


anticommute 


(A8) 


where we used anticommutation of THa and THb- We recover the result Eq. (A5) for NqAi = 2. 

In fact, one can view VH on a single orbital as fractionalization (projective representation) of the complex-conjugation operator 
K,. While /C is a bosonic operator squaring to 1, the same does not hold true in general for 7^. For iVorb = 1 (or 3), VH is a 
fermionic operator. For iVorb = 2, VH^ = — 1 and the system has a Kramers degeneracy. Only for iVorb multiple of 4 the VH 
operator behaves like complex-conjugation. Equivalently, we can only find eigenstates (with definite fermion parity) of the TH 
operator when iVorb S 4Z. When iVorb ^ 4Z and the Hamiltonian respects particle-hole symmetry, there will always be a ground 
state degeneracy. (Note that at odd iVorb all states break PH symmetry, so there is trivial two-fold degeneracy if the chemical 
potential is chosen such that the Hamiltonian is PH-symmetric.) 

The situation parallels the “BDI-wire” studied by Fidkowski and Kitaev,"^® where they considered a fermionic system which 
possesses a time-reversal symmetry squaring to -1-1. There is a Zg classification, i.e., phases described by an integer n defined 
modulo 8. When n ^ 8, the time-reversal operator also fractionalizes at each end of the ID system.^®’^* The case with VH oper¬ 
ator acting on a single orbital is akin to the Fidkowski-Kitaev n = 2 case. The T’H operator acting on two orbitals corresponds 
to the n = 4 case. When acting on four orbitals, the T’H operator squares to -fl, which is analogous to the n = 8 = 0 (mod 8) 
case in the BDI-wire. In the latter case, the wire is in a trivial topological phase and has no ground state degeneracy. 


Appendix B: Quantum Hall on a cylinder 


To make the manuscript more self-contained, here we describe our setup for the DMRG studies of quantum Hall problems 
on the infinite cylinder. We also describe symmetries in this setup, which are useful, e.g., in the discussion of putative effective 
field theories for the CFL 3.ti/ = 1/2. The electron kinetic energy Hamiltonian in the Landau gauge (A^, Ay) = (0, Bx) is 


Ha. 


kin. 


2m 


Vi 



(Bl) 


where (.b = \/fi./(ei3) is the magnetic length. Orbitals in the lowest Landau level in this gauge are 


Tj{^) 



(B2) 


where r = {x, y). The orbitals are labeled by an integer /; such an orbital is a plane wave in the y-direction with wavevector 
ky = {2tv/L y)j and is localized in the a;-direction around position Xj = £\(2'k / Ly)j . We can naturally view these orbitals as 
forming a ID chain labeled by j G Z. Going to the second-quantized language, we expand the continuum electron annihilation 
operator c(r) in terms of the annihilation operators cj for these orbitals, c(r) = ipj(r)cj, cf. Eq. (A2). 

On the cylinder, the above electron kinetic energy has the following symmetries, which we readily translate to transformations 
of the Cj fields, i/ei.kin. is invariant under translation by an arbitrary amount Ay in the y-direction, which becomes multiplication 
by a j-dependent phase factor for the cj fermions; 


A generator for this symmetry is 


T^[Ay] : cj . 


J-^1/ 


(B3) 


(B4) 


which is proportional to the “center of mass” position of the system of electrons on the ID chain of orbitals. 

^^ei.kin. IS also invariant under translation by a discrete amount Ax = i'^{2Tr/Ly) in the x-direction, accompanied by a 
gauge transformation that respects the periodic boundary conditions in the y-direction. This discrete step is precisely such as 
to accommodate one flux quantum in the swept area of LyAx, and is also the spacing between centers of neighboring orbitals 
described above. Naturally, this symmetry becomes translation of the ID chain of orbitals, 

T^[Ax = e%{2Tr /Ly)] : Cj -)■ Cj+i . 


(B5) 
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The physics in an infinite 2D system is invariant under spatial rotations. However, when restricted to the cylinder, we only 
have 180 degree rotation symmetry left, which is the same as spatial inversion r —> —r. In the above orbital language, this 
becomes 


/ : Cj —>■ c_j . (B 6 ) 

In some situations it is convenient to consider inversion about a midpoint between two neighboring orbitals, which can be also 
viewed as a combination of the above inversion and translation by one orbital, 

I' : Cj —^ C-i+i ■ (B7) 

Finally, iFei.kin. is invariant under mirror reflection in the a;-axis, {x^y) {x,—y), combined with time-reversal T for 

spinless electrons (i.e., complex conjugation). We denote this symmetry as M^T, which becomes in the orbital language 


MxT : Cj Cj , i ^ —i . (B 8 ) 

Thus, this symmetry is simply complex conjugation in the ID chain, which ensures that the electron Hamiltonian including all 
interactions is real-valued in this basis. DMRG can then use real numbers only, which simplifies numerical calculations and 
automatically maintains this symmetry. 

We assume that all electron interactions also respect the above symmetries. In the standard procedure, we project the Coulomb 
interaction into the lowest Landau level and obtain 




E 


E 

j n>0,m>0 


WmnCjCj+nCj+m+nCjj_^j_2n + H.C. 


(B9) 


with calculable Wmn- The and Ty translation symmetries are manifestly present as chain translation symmetry and center- 
of-mass conservation respectively; the inversion I is also already imposed in the above form, while the anti-unitary mirror M^T 
requires Wmn G R- These details of the setup are all implemented in the DMRG but are not important for the rest of this 
Appendix. The electronic Hamiltonian becomes more complicated when we allow tunneling into a second quantum Hall well 
or allow Landau level mixing, see Refs. 26 for details, but the symmetry analysis remains. 

Let us now consider particle-hole transformation in the lowest Landau level. It is expressed in terms of orbitals in Eq. (A3) 
and is a symmetry of Lfei.int. 3tv = 1/2. In fact, it is easy to see that any four-electron term WjkimCjclciCm + H-C; with distinct 
j, k, I, m, is invariant under PH. Note that the interaction amplitudes Wjkim can be arbitrary complex numbers and do not need 
to respect any of the spatial symmetries discussed earlier, as long as j, k, I, m are distinct. If some indices coincide, e.g., j = I, 
a PH-symmetric form would have cjcj replaced with CjCj — which can be interpreted as requiring specific relations between 
interactions and (possibly site-dependent) chemical potentials. On the other hand, requiring A/ = A/,jb/2 in a translationally 
invariant system automatically satisfies this. Note also that any six-fermion interaction will violate PH symmetry, so generically 
it is not a symmetry of a truly microscopic electronic Hamiltonian that includes physical effects such as tunneling to another 
layer or Landau level mixing. 


Appendix C: Numerical methods 
1. Infinite DMRG 

Our implementation of infinite quantum Hall DMRG is described in Ref. 26. Here we provide some additional details specific 
to simulating the = 5 state via DMRG. 

DMRG is a ID method. Using the cylinder basis described in App. B, we map the problem to a ID fermion chain with a 
basis labeled by the orbital occupations rij. Each DMRG run starts with some orbital product state (e.g., | • • • 1010 • • •)) and then 
optimizes the state to reduce the variational energy. As discussed above, the anti-unitary mirror symmetry M^T is maintained 
automatically by working with real-valued arithmetics. Furthermore, the DMRG explicitly preserves Ky, so the initial state 
fixes a sector with definite Ky. When started in one sector, the DMRG then resides solely in that sector. We need to check 
multiple sectors to determine which one contains the ground state. In the present study, we use a four-site unit cell when we set 
up MPS for the infinite-DMRG, and we can readily check that “root” configurations 1010, 0101, 1100, and 0110 generate four 
different such sectors. They also exhaust all different sectors accessible with such a unit cell at half-filling, since 0011 can be 
connected by center-of-mass-preserving terms to 1100 (e.g., on a torus with A),ib multiple of four), while 1001 can be connected 
to 0110. Note also that the 1010 and 0101 sectors are related by T^, and so are the 1100 and 0110 sectors, which gives exact 
two-fold degeneracy of eigenenergies. Such a degeneracy in the half-filled Landau level is a consequence of the commutation 
relation T^Ty = —TyT^ s.lv = 1/2. To be precise, the last equation is for an x Ly torus, where we consider translations 
Tx = /Ly)], which is the same as in Eq. (B5), and Ty = Ty[t\(2T: / Lx)]. Using the generator Ky, Eq. (B4), the latter 
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can be written as Ty = exp[i(27r/iVorb) where TVorb = L^Ly /{2'Ki’^g) is the total number of orbitals, which is the same 

as the number of sites in the ID chain representation. In our numerical calculations, we then restrict to studying just two sectors, 
0101 and 0110. 

Note also that the VH transformation Eq. (A3) and Ty anticommute (assuming even iVorb). Hence, genetically the VH will 
connect sectors with different ATj,; however, there are exceptions where the VH acts within a sector. For our four-site unit cell, 
the VH connects the 1010 and 0101 sectors. On the other hand, the VH maps 0110 to 1001 which belongs to the same Ky sector, 
so it acts within this sector; the PH is now a discrete symmetry for the Hamiltonian restricted to this sector. 

The DMRG can spontaneously break the PH symmetry due to quantitative errors induced in the results by the finite bond 
dimension. As bond dimension is increased, the energy approaches the infinite-bond-dimension limit from above as a power 
law in the bond dimension.^® The DMRG is trying to find the lowest energy state of the system, and it can choose between a 
“cat state” which is a superposition of the HER and HER states, or a symmetry-broken state which is only HER (or HER). The 
exact energy splitting would be zero or exponentially small in the cylinder length. However, for a finite bond dimension, the 
cat state must use half of the bond dimension to describe the HER and half to describe the HER. Therefore each piece has its 
bond dimension effectively halved, which introduces an algebraic increase in its energy. This energy increase makes the cat state 
energetically less favorable than a pure state. So if our DMRG-optimized state at finite bond dimension does not break the PH 
symmetry, we conclude that the true state in the infinite-length quasi-ID system does not break this symmetry. Note that this 
argument is actually valid for any discrete symmetry breaking in quasi-lD. 

Once we have the optimized DMRG state (matrix product state) at a given bond dimension x, we can calculate the correlation 
length ^ in this state from the subdominant eigenvalue of the transfer matrix. In a gapped system, ^ converges to the true 
correlation length in the ground state. In a critical system, it diverges with increasing x, while at a fixed x it can be viewed as 
an effective cutoff length set by that bond dimension. This is the length scale plotted in Fig. 3 vs entanglement entropy between 
left and right parts of the system, which is readily measured in the same MPS state. 

We also study an overlap of the optimized DMRG state with its particle-hole conjugate. Using similar transfer matrix tech¬ 
nique, we can naturally calculate an effective overlap per one orbital, 1 — e, used in the main text. 


2. Correlation functions 

Importantly for this work, we also measure density-density correlation function. We first specify our correlation function in 
general. We consider a fixed number of electrons moving in continuum in a region of volume Vol = L^Ly and define 

£)(r - r') = (: <5p(r)(5p(r'):) = (: p{v)p{r'):) - p^ . (Cl) 

Here &p{v) = p(r) — p is the deviation of the density from its average value p = A"e/Vol, and we are also assuming translational 
invariance. The corresponding structure factor is 


D{q) = JdrD{r-r')e = {:Spq6p^q:) = {6pqSp-q) - p , (C2) 

where Spq = J dr i5p(r)e“®‘^ '’/V^Vol. The last expression follows from simple manipulations (e.g., in the first-quantized lan¬ 
guage), and one can use it to argue that thus defined id(q) is a continuous function with limiting behaviors D{q —> c») = 0 and 
id(q —>■ 0) = —p. On the other hand, for numerical evaluations it is easier to use the original expression, which gives 

^(q) = (:PqP-q:)-p"Voldq,o , (C3) 

where Pq = fdr p(r)e“*‘i '’/VVol. (Note in this rather unusual convention, (po) = pVVol.) 

We now specialize to the quantum Hall problem on the cylinder. From the outset, we will assume that is very large and 
treat the corresponding wavevectors as continuous, while Ly is finite and the wavevectors qy are discrete. Using the basis of 
orbitals in Eq. (B2), the density operator becomes 


p(r) = c^(r)c(r) = 


^ V^U'-j)v g- 
3,3' 


(ri-X p^ + Cx-x-,)^ 






(x-X 


—'Hl’U 

V e y e 






La, 




(C4) 
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In the last line, Xj^rnl 2 = ^B(27r/Ly)(j + m/2) is a position half-way between orbitals j and j' = j + m. From the above 
expression, we immediately see that operator Om{j) = c'jCj+m carries transverse momentum qy = (27r/Ly)m. We can now 
calculate the Fourier transform of the density operator. 


27r \ e 

P \ <lx,qy = -j-m ) = 

L/t, 








Since p_q = we have 


PqP-q = 


_ 

e 2 

Vol 




-Jt,,) 




_ ^gq 
e 2 

Vol 

e 2 
Vol 




J,,? 


^orb E ^0'=™ (4cA+m)^ ^ 

A 


(C5) 


(C6) 


In the last equation, we assumed that the state is invariant under translations of the chain, and iVorb is the total number of sites in 
the chain, which is the same as the number of fluxes through the x Ly region. 

Now our correlation function becomes 


, 1 ^ '1 ( 27r ^ Vorb 









q,0 


iV^ 

e 

N, 


orb 


(C7) 


We have TVorb/Vol = l/(27r£^). Also, V//Vorb = NeV, and we can verify that at {q^, qy) = 0 the term in the square brackets 
is equal to —i/ and is actually independent of the number of sites in the chain (i.e., cylinder length). While not obvious from the 
last expression, our earlier discussion after Eq. (C2) implies that thus defined structure factor is a continuous size-independent 
function of q^ near q^ = 0 at qy = 0. We have D{qx —>■ 0, gj, = 0) = —iyl{27ri'^) = —l/(47r) in units used in the main text, 
and the DMRG measurements in Fig. 1 indeed give this value. 

The above expression can then be used in the infinite cylinder DMRG setup. The sum on A in principle runs from —oo to oo, 
but in practice summing over a few hundred orbitals provides sufficient momentum resolution. When plotting the results, we 
often omit the Gaussian factor, i.e., we show D{q} = 

Finally, we note that when we work with a given unit cell and fixed root configuration, the state is guaranteed to be transla- 
tionally invariant only under translations by unit cells. The sites inside the unit cell need not be equivalent. In the limit of very 
large Ly these sites are spaced by a very small amount £^(27r/Lj,), and such microscopic variation is irrelevant in the 2D limit. 
In our finite Ly cylinders, we simply average the position of the “0” site in Eq. (C7) over sites in the unit cell. 


Appendix D: Additional numerical data 

Here we present a number of additional numerical results which support the conclusions drawn in the main text. 

In Eig. 6 we map out the Eermi surface using the lengths Qm of the various “wires.” The ky coordinate of each point is given 
by 21 : 7 x 1 / Ly, with m integer or half-integer depending on boundary conditions, and the coordinate is determined from the 
singular wavevectors qx = Qm in the qy = 0 structure factors (see Eig. 2), namely kx = Qm/‘^- We see that the resulting points 
lie very close to a Eermi surface with kp — 1, with slight deviations which are needed to satisfy Luttinger’s theorem in the 
quasi-ID system. 

In Eig. 7 we show another view of the Ly = 24 density structure factor D{q) shown in Eig. 1. Here each line has a different 
qy, with a vertical offset proportional to qy. We can see both the accumulation circle at 2kF, and also other finer features 
representing other singularities present in the quasi-ID system. Note that for this wide cylinder, our largest bond dimension used 
still effectively cuts off the correlations at a length scale of order 10, see Eig. 3, which explains some rounding of the features 
compared to cylinders with smaller Ly. Note also that in this plot, as well as for all other sizes [cf. Eigs. 2(b) and 9], we do not 
see any features at = 0 and small qx. This is consistent with a much weaker singularity ~ \qx\^ expected in this case (see 
App. 13), which has continuous first and second derivatives. 

In the main text we measured correlators of the PH-even observable P(r). A natural observable to use for P{v) would be 
i5p(r)V^p(r). However, in the orbital basis each p contains two summations over many sites on the chain, so the full correlator 
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FIG. 6. Map of the Fermi surface by collating data from various circumferences [same systems and color scheme as in Figs. 2(c) and 3]. 
Note that these points deviate slightly from the 2D circle, in order to satisfy Luttinger’s theorem for the quasi-lD system with discrete wires, 
cf. Fig. 2(c). 
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FIG. 7. Alternative 3D view of the data in Fig. 1 for the Ly = 24 cylinder with eight-wire CFL. We show density structure factor measured 
at fixed qy ranging from —8(27r/Ly) to 8{2n/Ly) from bottom to top (labeled on the right axis in units of 2njLy}. Each curve has a vertical 
offset proportional to Qy. While we see the appearance of the two-dimensional 2k f accumulation circle from the low-energy composite 
fermion excitations, we also see finer features corresponding to singularities in the quasi-lD system with finite Ly. 


contains thousands of terms. Instead of calculating all of these, we calculated only a few dozen such terms, reasoning that all 
terms should contain the same singularities. None of the terms we calculated shows a singularity at the exact backscattering 2kF. 
In Fig. 5 we plot (P 5 (q)P 4 (—q)), where P 5 = cjcj+inj+s, P 4 = njCj+ 3 cj_|_ 4 . Note that P 5 carries qy = 2 t:/L y and P 4 carries 
the opposite qy, and the exact backscattering probed here is between the right-mover composite fermion at ky = jLy and 
left-mover at ky = —^/Ly, see bottom picture in Fig. 2(a). We chose this correlator because it clearly shows the reemergence of 
the 2kF exact backscattering peak when the particle-hole symmetry is broken by hand, but other correlators show this behavior 
as well. In Fig. 8 we take a detailed look at this correlator in the system with PH symmetry and identify all singularities with 
various processes, using information about the lengths of the wires determined from Fig. 2(b) in the main text. We have labeled 
many singularities, including those for all allowed two- and four-body processes transferring composite fermions from the left 
to right side of the Fermi surface (with total qy transfer of 27r(Ly), but we do not see the exact backscattering 2k f singularity, 
as expected in the presence of the PH symmetry. To show that our results do not depend on the particular PH-even observables 
presented, we also show the correlator (P 2 (q )^2 (~q))’ where P 2 = c]cj+inj-|_ 2 . 

Most of the Fermi surface evidence in the main text was shown for Ly = 13 realizing four-wire Dirac-CFL, since here the 
finite-entanglement errors are the smallest. However, we have similar evidence for several other sizes. In Figs. 9 and 10 we 
show the same data as in Figs. 2(b) and 8 , but for Ly = 16 realizing five-wire Dirac-CFL. The singularities are more rounded at 
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FIG. 8. Derivative of correlation functions for PH-even operators as in Fig. 5 but for a single layer. The cylinder has Ly = 13 and realizes 
four-wire CFL (APBC; 0110 root configuration). The two curves correspond to different correlators used, (P 5 (q)P 4 (—q)) (same as in Fig. 5) 
and {P 2 (q)P 2 (—q)), see text for details. We measure at qy = 27r/Lj, since it can probe exact composite fermion backscattering with such 
boundary conditions. We can label all features by transfer processes contributing to such an operator P(r). We see one-fermion transfer 
A —>■ B', while the exact backscattering P —>■ G' is clearly missing, as expected for the Dirac CFL. Amazingly, we also see higher-order 
processes corresponding to transferring two left-moving composite fermions to the right-moving points, which we labeled AD —> B'D', 
BD —>■ C'D', AC —>■ B'C' (in fact, these exhaust all distinct wavevectors for such four-fermion terms with qy = 2TT/Ly). We even see a 
contribution from a six-fermion term transfening ABC —>■ A'B'D'\ the only other six-fermion term that transfers three left movers to three 
right movers would be ABD —>■ A'C'D', which however is a combination of three exact backscatterings and hence is odd under PH and 
cannot contribute to PH-even observables. 



Momentum 

FIG. 9. Same as Fig. 2(b) but for Ly = 16, which realizes five-wire CFL as shown (PBC for the composite fermions; root configuration 1010). 
Note that the singularities are more rounded here since the largest cutoff length ^ accessible with our finite MPS bond dimensions is about 
four times smaller than in the Ly = 13 case, cf. Fig 3. The top panel shows fractional derivatives whose powers are chosen individually for 
each singularity, roughly to turn it into a step singularity. We estimate the lengths of the wires from the qy = 0 data, and the dashed lines are 
determined from those estimates. 
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FIG. 10. Same as Fig. 8 but for Ly — 16 and transverse wavevector qy — 0, plotting derivative of the correlation function of operator 
Pi = hjhj+i. This cylinder is the same as in Fig. 9 and realizes five wires as shown. PH-even observables cannot obtain contributions 
from precise 2kp backscattering, which here corresponds to C ^ C' transfer of one composite fermion. On the other hand, A ^ A! and 
B ^ B' are not exact backscatterings in the 2D sense, so these processes do not have definite PH quantum number, and we do see features 
corresponding to these transfers. We see many other features, e.g., from four-fermion terms involving transferring two left movers to two right 
movers, and even from six and eight-fermion terms. In fact, we see all possible distinct wavevectors from such four-fermion processes and 
mark them on the plot, and we also mark the most visible higher-order processes. 


Ly = 16 compared to Ly = 13 due to more significant finite entanglement effects, i.e. the effective cutoff length ^ imposed by 
our finite MPS bond dimension is smaller at these sizes, see Fig. 3. 

In Fig. 10, the PH-even operator whose correlator we measure is simply Pi = Tijrij+i, and because of the PBC for the 
composite fermions, we probe the exact backscattering by measuring at qy = 0. We clearly see singularities at wavevectors 
corresponding to processes transferring one composite fermion from A to A! and from B to B', which do not have dehnite PH 
quantum number and can hence contribute to both PH-even and PH-odd obserwables. On the other hand, there is no singularity 
corTesponding to transfer from C to C, in agreement with the absence of exact backscattering. This figure also shows that our 
measurements can pick up signatures coming from contributions from higher-order processes. In fact, we can identify all six 
distinct q^. wavevectors corresponding to processes transferring two left-moving composite fermions to the right, with the total 
qy = 0. Each such distinct wavevector is labeled by one process with such q^, e.g., AB A'B', while there can be multiple 
processes with the same {qx,qy = 0). We even see wavevectors corresponding to processes transferring three and four left- 
movers to the right, but we have not tried to identify all of them systematically. [Note that six-fermion terms ACE A'C E' 
and BCD B'C'D' are again combinations of three exact backscatterings and are not allowed in the PH-even operators.] It is 
amazing how much information the DMRG gives us, and all is consistent with the proposed quasi-ID descendants of the Dirac 
CFL! 

In Fig. 11 we show how the overlap behaves as a function of bond dimension, using a modified Coulomb interaction (1 — 
x)Vo + xVi- Here Vq and Vi are the (Gaussian-cutoff) Coulomb interaction projected into the lowest and Landau level, 
respectively. At a; = 0 the system is in the CFL phase, while at a; = 1 the system is in the iz = 5/2 state with spontaneous 
symmetry breaking.'^ We see that in the CLL phase the overlap is over one hundred times smaller than in the symmetry 
broken phase, and that it decays to zero as a power law. 


Appendix E: The effect of Landau level mixing on 2kF backscattering 

Several experimental perturbations break PH symmetry, including disorder, incomplete polarization of the electron spin, and 
Landau level mixing. Here we comment on the last two effects. 

At many of the gapped plateaus, a combination of a Zeeman energy and interaction-driven ferromagnetism fully polarize the 
electrons. However, at ^ = 1/2, the polarization in fact depends sensitively on the electron density and the strength of the 
H-field used to achieve half-hlling. Lor spin-unpolarized electrons, LLL particle-hole transformation exchanges zz o 2 — zz, as 
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FIG. 11. A plot of e vs. bond dimension, where the overlap (per orbital) between our wavefunction and its particle-hole conjugate is 1 — e. 
Data was taken at Ly = 16 using a modified Coulomb interaction {1 — x)Vo -\- xVi, where Vo is the Coulomb interaction in the lowest Landau 
level and Vi is the interaction in the L* Landau level. When x is small the system is in a CFL phase with particle-hole symmetry, while when 
a; ~ 1 the system is in the Moore-Read phase and spontaneously breaks this symmetry. We see that that at the largest bond dimension used, 
when X is small e is two orders of magnitude lower than when the symmetry is broken, and it decays as a power law as the bond dimension is 
increased. The values quoted in the main text are obtained by fitting these curves to a power law. 



FIG. 12. Effect of Landau level mixing on the 2kp backscattering. Shown here is the derivative of the correlator (P 5 (q)P 4 (—q)) at Ly = 13. 
We see that including an additional Landau level at mixing parameter k does not lead to a reappearance of the 2k f singularity even for 
unrealistically large k (see text for details). The main plot shows data at k = 0 and k = 3, and we see that the two plots have the same 
features. The inset focuses near the backscattering wavevector and shows data up to the value ac = 5, above which there is a transition to a 
gapped phase. 


it must act on both spin species, so the particle-hole symmetry discussed in the main text is broken. In the CF picture, the CFs 
are spinful, and will polarize depending on the ratio of the CF Fermi energy to Zeeman energy, Ep/Ez- Note that Ez oc B, 
while Ep ~ e^/(47re£B) oc cfB, since the CF Fermi surface is due entirely to Coulomb repulsion. Thus at high fields the 
Zeeman energy dominates and the CFs fully polarize. This transition is accessible experimentally^^ thus in order to preserve 
PH symmetry, we assume the high-field side. 

Particle-hole symmetry is exact only when the Hamiltonian is projected into the LLL, which is justified when the cyclotron 
gap Ey; = fujJc oc P is large. However, the Coulomb scale Ec = ^ j{E kcIp) is not negligible, so the LLL projection is violated 
by terms of order k = Eq/E^ oc with k ^ 0.4 in typical LLL experiments. 
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We can ask whether the 2kF backscattering returns for experimentally relevant k, as it did in the bilayer setup in the main text. 
To investigate, we project the Coulomb Hamiltonian into the two lowest Landau levels, n = 0 and 1, and tune k by changing the 
cyclotron energy. The result is shown in Fig. 12. At k = 0, the measurement is identical to the ^ = oo curve of the main text. 
Fig. 5. 

The singularities away from 2kp change quantitatively with k, but the 2kF backscattering does not reappear, in contrast to 
the bilayer numerical experiment in the main text. Presumably there is some 2kF singularity with an amplitude too small for 
us to resolve. Thus, according to this operational definition, PH-breaking is “weak.” Furthermore, numerically we find that at 
K = 1 the occupation per flux in the n = 1 LL is fi « 0.003. We observe such “weakness” of the effective PH-breaking over 
the entire range of the gapless CFL phase in this numerical experiment, including unrealistically larger k, while for even larger 
K > 5 our model goes into a gapped phase (likely the Moore-Read phase). 


Appendix F: The Dirac CFL 

Here we briefly review the proposed Dirac-CFL theory in Ref. 20. Like in the TI, the composite fermion '0 cf;s is a two- 
component fermion (s = 1, 2) on which PH acts as = icr^i^ce- The component “s” is not the electron spin, which is 

polarized, but is rather an emergent degree of freedom which can be pictured as the dipole moment of the CF. The new ingredient 
in the Dirac-CFL, compared to the surface of a TI, is the emergent gauge held a^. The = 1/2 phase is then proposed to be 
equivalent (dual) to (2 h- 1)D quantum electrodynamics: 

>Cqed3 = + • ■ • (FI) 

Here '0 = ' 0 ^ 7 °, with 7 ° = cr^, 7 ^ 7 ^ = cr^, 7 ^ 7 ^ = cr^; we have also set velocity of the Dirac fermions to 

1 and used common writing of the coupling constant in the Maxwell term /^. In addition to the Dirac and Maxwell Lagrangians, 
the last term “dopes” the CFs away from the Dirac point to finite density, as it is equivalent to a static background charge B/Af. 
Thus, like in HLR, the CFs have a Fermi surface of radius kp = and are also coupled to a dynamical gauge held. 

To probe the physical electron, we introduce the physical vector potential A. It appears in Eq. (FI) via the replacement 
^ + ^)dyA\. The following relation for the electron density follows: rie = ^V x a -f 

The Dirac CFL theory maintains the particle-hole symmetry and thus differs from the original HLR construction. However, 
some crude properties are quite similar in the two theories. For example, both theories predict that in the case of Coulomb 
interactions, the leading non-analyticity in the long-wavelength density structure factor is |qp log(l/|q|). While this is a very 
weak singularity and difficult to see in numerical simulations, it would be interesting to attempt different scaling analyses of the 
DMRG data to study such non-Fermi-liquid aspects of the CFL in 2D. 


Appendix G: Locking of momentnm with CF boundary conditions 


We observe that the BC of the CF is locked to the momentum Ky = jTij. While the absolute value of Ky 

is difficult to define on an infinite cylinder, we can safely consider the momentum per electron relative to (say) the orbital 
configuration • • • 0101 • • •. Our observation of PBC for 0101 and APBC for 0110 is consistent with the relation 


AiT, 


= N, 


2 f T>int 

Ly 27r ’ 


(Gl) 


since the difference between the fl)int = 0 or tt orbital configurations is such that every second electron is bumped by one orbital, 
with AKy = . We now derive this relation within the Dirac-CFL theory. 

Note that any discussion of momentum at finite held is subtle as translations obey a magnetic algebra: 


T,ix)Tyiy)T,{-x)Ty{-y) = . (G2) 

On a cylinder of circumference Ly, translations along the length of the cylinder are broken down to a discrete subgroup T^{x = 
On a torus, translation is also broken along y, and we have 


T, 





(G3) 


Thus when the dual theory, finite density QED 3 , is placed on a cylinder / torus, we expect that translation symmetry must be 
broken in some subtle fashion. We will see this is indeed the case. 
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Let us begin with the dual theory on a torus in the presence of a uniform background magnetic field B. At the end of the 
calculation we can take a limit to recover the infinite cylinder result. We also take the temporal direction to be compact. The 
Lagrangian is 


L = (G4) 

Note that there are several ways to rewrite the CS term (related by integration by parts) and various boundary conditions one 
may place on and 4}. Let us stick to the form of the Lagrangian in Eq. (G4). 

We discuss the physics in a Hamiltonian formalism. We consider the situation when the physical Landau level is half-filled, so 
there is no net flux of {d^ay — dya^) through the spatial torus. We take all fields (a^, ip) to be periodic along x and y directions. 
We have large gauge transformations, 


UxttxUl =ax + UxayUl = a^, Uxip{x, y)Ul = y) , 

UyaxUl=ax, UyayUl=ay + —, Uyip{x,y)Ul = e''^'^y/^'>ip{x,y) . 

Ly 

As discussed in Ref. 22, for a single Dirac cone regularized in a time-reversal invariant manner, 

UxUy = -UyUx . (G 6 ) 

This is the parity anomaly. We require the ground state wavefunction to carry Ux = ^ and (7^ = 1; equivalent results are 
produced if we exchange the role of x, y. 

Now, let us discuss the translational symmetry. The key issue has been previously identified as a feature of QED 2 (the 
Schwinger model) at finite density: translation symmetry is anomalous, and the 6 >-vacua are permuted by translations. 
When QED 3 is placed on a cylinder, an analogous effect occurs. When B = 0, we have the conserved gauge-invariant energy 
momentum tensor. 


= 1pi^^^D''1p 


1 

4e^ 


pti\pi^ _ 


- iax)ip 


1 

4e^ 



(G7) 


We have not symmetrized the energy momentum tensor as symmetrization will not be important for our purposes. When 5 = 0, 
we have dyT^'^ = 0. If one couples the theory to an external current, 

L = ipi'yf^{df, - iay)ip - + a^ierct (G 8 ) 

then the energy momentum tensor is no longer conserved: 

= r^jx,ext (G9) 


This result follows from classical equations of motion. It does not rely on the Dirac dispersion of the fermion (would be still 
true for a non-relativistic fermion, provided we use an appropriately modified Now, let’s consider a uniform background 

charge density j^xt = {—p, 0). In our case, p = ^. Then, 


= 0 , 


(GIO) 


Thus, energy is conserved, but momentum is no longer locally conserved. We may try to fix the problem by redefining. 


foz ^ pm p pij ^ p^j _ 


(Gil) 


We now have = 0. However, is now no longer gauge-invariant. The situation appears to be a bit better for the global 

momentum built out of 


pr 



J (fx (r°* -f pm) 


(G12) 
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The momentum P® is conserved, moreover, it is invariant under small gauge-transformations with a{x) periodic in 

space). However, it is not invariant under large gauge transformations (G5), 


U^P^Ul = P^ + 2 tt pLy = P^ + 




u^pyul = py, 


UyP^ul = Pp Uypyul = py + = py + 


'^^orh 


(G13) 


where A^oib = BLxLy/2'K is the number of flux-quanta. Thus, there are no infinitesimal translations which commute with the 
Hamiltonian and with the large-gauge transformations and Uy. However, there are discrete translations which do. Define 




(G14) 


Ty{by) commutes with Ux for any by. It commutes with Uy provided that 

by = my eZ 

Norh 


(G15) 


This is the correct quantization of translations in the Landau level, Eq. (G3). As for Txibx) it commutes with Uy for any bx, and 
with Ux when 


^ AT 

-iVorb 


mx € Z 


We know that in the Landau level, smaller translations by bx = should also be a symmetry. However, we have, 

f Ux = -UxTx ( 


Tx , 

v^orb; 

To get rid of the minus sign, consider a modified translation operator. 


V -^orb / 




(G16) 


(G17) 


(G18) 


We know that Uy is a symmetry of the Hamiltonian, so Tx is also a symmetry. Moreover, due to Eq. (G6), Tx now commutes 

since {7^ = 1 on the physical Hilbert 


with Ux (and with Uy). Note that on the physical Hilbert space T^ ~ ( 


space. 

Einally, let us compute the commutator of Tx 




(n 1) (wt.) 

P® = / (fxT^^ 


. Defining the bare, unconserved momenta as 


for any local operator 0{x) we have, 

[p®,o(*)] =d,o{x). 

Moreover, since the action contains no Chern-Simons term for a^,, Gx and ay commute. Therefore, 

pppy] = -ip f d^x{dxay - dyGx) = -IB{N, - 7Vo,b/2) 


(G19) 


(G20) 


(G21) 


where N^. is the number of physical electrons in the Landau level. Strictly speaking we are considering the case where there 
is no background magnetic flux of a (A7 = TVoib/2), so the commutator is 0. (However, further consideration indicates that 
Eq. (G21) continues to hold even away from half-filling). Then, 


Tx 



Norh 


= e 


27711/ 


T, 


Norb 


Tx 


Lx \ 

N,,b) 


(G22) 


as required by the magnetic algebra, Eq. (G3). In particular, at iz = 1/2, Tx and Ty anti-commute. 
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Now we are ready to understand the change in momentum when we have an “even” vs “odd” number of CF modes on the 
cylinder, i.e., when Oy = 0 vs Oj, = In this case, the expectation value of the un-corrected momentum is {py) = 0, since 
the Fermi sea fills in a rotationally symmetric fashion. Thus the corrected momentum is 

{py)=pj(fxay (G23) 

Hence, the difference in between Uy = 0 and ay = ^ is, 

A{py) = ^^ (G24) 

Ly Z 


in agreement with the hypothesis of Eq. (Gl). 

One may ask if this difference is meaningful on a torus (rather than a cylinder). Consider the operator Ty{2Ly/Norh)- This is 
the double of the minimal allowed translation operator on the torus. This operator acting on our two states gives results differing 
by a — sign. So, indeed, the difference in momenta is meaningful on a torus. Further note that since we are using the double 
of the minimal translation to distinguish between the two states here, the subtlety involved in the definition of Tx does not enter 
here (i.e., there is no difference between the two directions). 


Appendix H: Half-fllled Landau level and topological superconductor in class AIII 

In this appendix, we elaborate on the connection between the half-filled Landau level and the surface state of a 3D topological 
superconductor (TSc) in class AIII.^°’^^ TSc’s in class AIII in 3D are close cousins of the familiar 3D topological insulators. 
Like the TIs, they are phases protected by a combination of (7(1) symmetry and time-reversal T. However, unlike in the TI, here 
the (anti-unitary) time-reversal symmetry acts on the charged fermion i/) in a particle-hole manner, inverting the (7(1) charge of 
■0. One may physically think of the (7(1) symmetry as conservation of the ^-component of spin in a superconductor, thus the 
nomenclature. Unlike non-interacting 3D TIs, which have a Z 2 classification, non-interacting 3D TSc’s in class AIII have an 
integer classification. The surface of a phase in class 1 / G Z supports \v\ Dirac cones. We focus on the 1 / = 1 phase with a single 
surface Dirac cone. 


H = [—((9a; — ieAx)a^ — i{dy — ieAy)a'^^ 0 (HI) 

Here, 0 is a 2-component fermion, which transforms as T : 0 —> <7^0^, (7(1) : 0 —e*“0. For future convenience, we have 
included the coupling to a (7(1) gauge field A. 

Let us place the surface in a magnetic field B = d^Ay — dyA^. Unlike in a TI, here the magnetic field does not violate the 
time-reversal symmetry, due to the particle-hole nature of the latter. The surface spectrum breaks up into Landau levels with 
energies, En = sgn(n)'(;-\/2e|i?| |n|, n = 0, ±1, ±2,.. .. T-symmetry forces the n = 0 LL to be half-filled. If one takes the 
limit where the Landau levels with |n| > 1 decouple, then the surface theory becomes identical to the 2DEG in the LLL level. 
In particular, the T symmetry of the surface maps to the PH symmetry of the LLL. However, unlike in a 2DEG, here the entire 
spectrum is T-symmetric, i.e., there is no need to project into the n = 0 LL for the PH symmetry to emerge. Notably, the 
time-reversal symmetry acts on the full bulk Hilbert space of the 3D TSc in a local manner, unlike the non-local action of the PH 
symmetry in the lowest Landau level of a 2DEG. Conventional wisdom states that surfaces of 3D topological phases cannot be 
imitated in a purely 2D setting. However, this wisdom assumes that the 2D implementation has a local symmetry action, while 
the PH symmetry of the 2DEG in the LLL is non-local, as explained in the main text and in App. A, which allows it to mimic 
the T-symmetric surface of a class AIII TSc. 


Appendix I: Field theory for the quasi-ID descendant of Son’s Dirac CFL 
1. Fields near the Fermi surface for Dirac composite fermions at finite chemical potential 

In this Appendix, we use Son’s Dirac CFL theory^® (reviewed also in App. F) to propose an effective field theory for the 
quasi-lD descendant states on the infinite cylinder. In the same spirit as in Son’s paper, we postulate that low-energy degrees 
of freedom are some fermions (loosely referred to as “composite fermions”) coupled to a dynamical gauge field. Here we do 
not derive these fermion fields microscopically (but see ideas in Appendices H and J). Instead, we follow Son and postulate 
transformation properties of the fields under all microscopic symmetries. This constitutes in principle a complete specification 
of the theory and allows in particular identification of low-energy field contributions to any observable. 
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Son’s theory is formulated as a Dirac theory at a finite chemical potential. Here we will focus on fields residing near the Fermi 
surface where the Fermi level cuts through the Dirac cone. We encode the Dirac cone nature in the symmetry transformation 
properties of the fermions near the Fermi surface. We do this rather than keeping the two-component Dirac field, since this 
approach translates immediately to the quasi-ID setting. 

The Dirac Hamiltonian v{a^kx -f o'^ky) — fjL has eigenergies ±w|k| — p. Assuming /r > 0, the low-energy excitations come 

from the “H-” branch and reside near the Fermi surface at |k| = kp = fi/v. We take the corresponding eigenvector as ’ 

where ctk denotes an angle formed by the vector k with the fc^^-axis. This fixes our choice of the phases of the wavefunctions 
for different k, which in turn fixes the transformation properties of the low-energy fields residing near the Fermi surface (note in 
particular discussion of the electronic particle-hole symmetry below). Specifically, we can expand Son’s two-component Dirac 
fermion field 'kcF(r) in terms of the low-energy fields residing near the Fermi surface, 

^CF(r)- fgiik] e**" Vk(r) +... , (II) 

k near the Fermi surface ' ' 

where /k(r) is a slowly varying field describing fermions near a Fermi surface patch at k, and we have also omitted contribu¬ 
tions from high-energy fields. The electronic quantum Hall problem has spatial symmetries and the anti-unitary particle-hole 
symmetry (in the lowest Landau level) described in App. B. We have in mind the same cylinder geometry in the Landau gauge 
and focus only on symmetries that are present in this geometry. For all symmetries, we use Son’s postulated transformations of 
the Dirac field T'cF(r) to obtain transformation properties of the fields near the Fermi surface /k(r) defined above. 

In the quasi-ID setting where the cylinder is infinite in the x-direction but finite in the y-direction, the low-energy fermions 
reside near Fermi points where the “wires” at discrete ky cut through the 2D Fermi surface, cf. Fig. 2(a). Let j label these Fermi 
points; kj is the wavevector at j; and Vj = 'Ua;(kj) is the group velocity along the cylinder and can be positive or negative 
corresponding to right- or left-moving fermions. The kinetic energy of the fermions is 

Hf,kin. = X! J dxf]{x){-ivjdx)fjix) . (12) 

3 

As described in App. B, the electronic model on the cylinder in the chosen gauge is invariant under arbitrary translations Ay 
in the y-direction. It is also invariant under a discrete translation Ax = i'^{27r/Ly) in the x-direction supplemented by a gauge 
transformation. When we say that we work with a long-wavelength field carrying momentum k, we implicitly understand the 
following transformation properties under the above symmetries; 

Ty[Ay] : /k ^ , (13) 

T,[Ax] (14) 

The microscopic model has the anti-unitary mirror symmetry combining mirror {x, y) —>■ {x, —y) and time reversal, which 
we denoted M^T in App. B. Son postulates its action on the two-component Dirac field as M^T : 'I'cF(a:, y) o’^'I'cF(a^, —y) 
(note that he denotes this symmetry as “VT,” while we preserve notation from App. B). Translated to our long-wavelength fields 
near the 2D Fermi surface, it reads 


Mx'T : f(k^,ky) f{-k^,ky) ^ i A (15) 

In particular, this symmetry implies that the Fermi surface is invariant under reflections in the fc^-axis and Vx{—kx,ky) = 

Xx{kx 1 ky ). 

The physics in the 2D system is invariant under spatial rotations. However, as in App. B, on the cylinder we only have 180 
degree rotation symmetry left, which is the same as spatial inversion r —^ —r. It transforms the Dirac fermions as / : T'cF(r) —> 
cr^'kcF(~r)j hence for the long-wavelength fields near the Fermi surface 

/ : /k ^ /-k . (16) 

This symmetry implies that the Fermi surface is invariant under inversions in fc-space and Ux(—k) = —z; 2 ;(k). Combined 
with the MxT symmetry, we see that four Fermi points {±kx,iky) are symmetry-related, with group velocities satisfying 
Vx{kx,ky) = Vx{kx,—ky) = —Vx{—kx,ky) = —Vx{—kx,—ky). This justifies our implicit assumptions about the composite 
fermion Fermi surface in the main text, cf. Fig. 2(a). 


2. Particle-hole symmetry and absence of 2kF backscatterlng 

We now consider the anti-unitary particle-hole symmetry present in the lowest Landau level at i/ = 1/2. Son postulates its 
action on the Dirac composite fermions as TH : T'cF(r) ——*c’'^T'cF(r), which translates for the long-wavelength fields near 
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the Fermi surface as 


e*“V-k , i^-i. (17) 

(Note that Son denotes this symmetry as “CT”, while we preserve label “TH” used everywhere else in our paper.) Note in 
particular that since a_k = ak + n, applying this symmetry two times takes any fk to —/k, i.e., T’H squares to —1 when acting 
on an odd number of fermions. This connects with our discussion of formal aspects of particle-hole symmetry in the main text 
and App. A. 

Furthermore, we observe that 

m : Aflf^k ^ -{Aflf^k)^ (18) 

for any complex coefficient A. This has an immediate consequence that particle-hole-even observables cannot have contributions 
from such “2kp” fermion bilinears. In the context of the surface of the 3D TI, similar considerations of a Dirac cone protected 
by time reversal symmetry imply absence of exact back-scattering by time-reversal-invariant perturbations. While the DMRG 
setup allows us to measure only properties of the ground state, we can nevertheless detect such physics by studying correlation 
functions of particle-hole-even observables and noting the absence of the corresponding 2fci? features, which we loosely refer 
to as the absence of the exact back-scattering. As we will see shortly, such 2fci? bilinears are enhanced and are prominent 

operators in the CFL theory, so their absence in the particle-hole-even observables is a dramatic feature. (Note that we do not 

claim complete absence of any singularity at 2kF- One can construct particle-hole-even contributions carrying such momenta, 
but only using more fermion fields, e.g., by multiplying f^f-k by a particle-hole-odd operator carrying zero momentum such as 
/k'/k' ~ /lk'/-k'- Such higher-order contributions are of course less important than bilinears would be.) 


3. Bosonized quasi-lD theory, gauge fluctuations, Amperean-enhanced 2kp bilinears, and density structure factor at long 

wavelengths 

The postulated Fermi points and the above symmetries fix the quadratic part of the Hamiltonian to the form in Eq. (12) with 
the symmetries of the Fermi surface as described above. Anticipating eventual multi-mode Luttinger liquid description in the 
quasi-ID limit, we want to use hydrodynamic bosonization approach to represent the fermionic system, and this will also allow 
us to readily treat gauge fluctuations. For the discussion of the particle-hole symmetry, it is convenient to group Fermi points at 
k and —k into right- and left-moving fields. Let m label such a pair of fields, which form together one gapless mode; henceforth 
m labels such “modes” before including effects of gauge field fluctuations in the CFL theory. We denote the right-mover of this 
pair as j — Rm and left-mover as Lm, and bosonize 

fpm = , for P = R/L = +/- , (19) 

with canonically conjugate boson fields 6m'ix')] = i7r6m,m'Q{x — x'), where 0(a;) is the Heaviside step function. Here 

rjm are Klein factors taken to be Majorana fermions, {?7m, Vm'} = 2(5m,mS which ensure that all fermion fields anticommute. 
The slowly varying fermionic densities are simply 


Pf,Pm = flmfpm = d,{Pcj)^ + e^)/{ 27 r) . (110) 

The kinetic energy in Eq. (12) now has sum over cumulative index j = Pm, with vpm = —VLm = Vm (and we assume Vm > 0). 
[Note that, alternatively, we could have grouped points {k^, ky) and {—k^, ky) into a right- and left-moving pair, since these also 
have exactly opposite group velocities along the cylinder. The physics does not depend on such choices, but the discussion of 
the particle-hole symmetry is slightly simpler with the made choice grouping k and —k used from here on.] 

Working in the imaginary time path integral, a bosonized Lagrangian corresponding to such iT/,kin in Eq. (12) is 

E/,kin = y] 7^ ■ (HI) 

^' Ztt tt 

m 

In the CEL theory (either HER or Dirac-CEL), the fermions are coupled to a dynamical gauge held, whose net effect in the 
quasi-lD system is to pin the overall “gauge charge” mode, thus reducing the total number of gapless modes by one.^^’^'^’^^ 
Indeed, working in the (1 h- 1)D space-time continuum theory in the gauge with only temporal component of the gauge held, 
the Lagrangian for the gauge held has an “energy” term proportional to {dxttrY and is coupled to the fermions via iUrPf,tot, 
where p/,tot = J 2 pm Pf,Pm = Thm ^x^ral'P- Integrating out the held produces a mass term for the overall gauge charge 
mode Qjxot = where is the number of R-L pairs (same as the number of wires). The remaining gapless 

modes can be obtained as linear combinations of that are orthogonal to Of ^ot, but we will not need details of these for the 
general discussion below. 
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One consequence of the above physics is that fermion bilinears that involve transfer from right-movers to left-movers are 
enhanced compared to the “free-fermion” theory. Indeed, consider such a bilinear 


ft 

B.m 


fhr. 


Vm Vn 




( 112 ) 


The combination of the 9 fields in the exponent has a component onto Of, tot, and since the latter is pinned, the fluctuating content 
in the exponent is reduced. This is quasi-ID manifestation of so-called Amperean enhancement of the 2k p bilinears in the 2D 
gauge theory: CF particle and hole from the opposite sides of the Fermi surface have parallel gauge currents which experience 
Amperean attraction. It explains our DMRG observations in the main text and App. D where bilinears (and also quartic terms) 
involving such transfers from right-movers to left-movers are prominent in the structure factors. 

Finally, we consider the behavior of the density structure factor at qy — 0 and small q^. Consider an operator in the long- 
wavelength theory dx[J2j where kj,y is the j/-component of the Fermi momentum at Fermi point j. This 

operator can be thought of as a derivative of a local momentum in the j/-direction (i.e., local piece of the conserved Ky). Since 
according to Eq. (B4), Ky also gives the ^-component of the total electron polarization operator on the cylinder, we see that 
the above operator gives effectively the electron density at gj, = 0 (i.e., averaged over the circumference of the cylinder). [We 
remark that while it is easy to verify that the proposed operator indeed has correct transformation properties under I, VH, and 
MxT, it is important to have the above connection to a microscopically conserved current expressed in the long-wavelength 
theory. Note also that unlike the 2D, there is no gapless gauge field left in the quasi-ID system, and the above contribution 
from the CF fields appears to be the dominant one.] We can now calculate contribution to the density structure factor and obtain 
D{qx, qy = 0) ^ |ga:P at Small q^. This is a very weak singularity, where only the third derivative becomes discontinuous. The 
density structure factors measured in the DMRG do not show any features at gj, = 0 and small q^, cf. Figs. 1, 2(b), 7, and 9, 
which is consistent with such weakness of the singularity. 


4. Stability of the Dirac CFL phase and irrelevance of partlcle-hole-breaklng perturbations 

Besides the fermion kinetic energy, the full theory also includes quartic fermion interactions (terms with more fermion fields 
or with derivatives are expected to be less important). The momentum-conserving quartic terms can be divided into two groups. 
The first group is forward scattering interactions 

^k,k'./k/k/k'/k' ^ Fpm,P'm' Pf,PmPf, P'm' 7 (113) 

with real-valued Fpm,p'm' as required by the Hermiticity. Under the anti-unitary mirror symmetry Pf,pm Pf-p-m (using 
convention where Fermi point —P, —m has the same ky as Pm but opposite k^), while under the spatial inversion Pf,pm —> 
Pf-p,m- We can easily write down corresponding conditions on Fpm,p'm'- The particle-hole acts on Pf,pm identically to 
the inversion, and since Fpm,p’m' are real-valued, it does not introduce new conditions. In any case, the forward scattering 
interactions are strictly marginal and only modify Luttinger parameters but cannot gap out the modes. 

The second group is Cooper (or backscattering) chanel interactions, 

Vk.k'/^/lk/-k'/k' -- Vm,m'flmflmfLm'fRm' , (114) 

with in general complex-valued The Hermiticity imposes Vm',m = Fm m" The anti-unitary mirror symmetry imposes 

V-ra-m' = m'’ while the inversion symmetry is automatically satisfied for any On the other hand, the anti-unitary 

particle-hole symmetry imposes additional conditions Vm,m' = where k = kpm is the wavevector of the 

right-mover in our m-th mode, and similarly for k'. Thus, the particle-hole symmetry essentially fixes the complex phase of 
Vm,m' to be 

Focusing on the CFL phase, we can make several interesting predictions. First, we note that the Cooper channel interactions 
can be all irrelevant. For example, in the theory with only the bare £/.kin supplemented by an infinitely strong gapping out 
condition on the overall gauge charge mode, i.e.. Of,tot = const, and no forward scattering interactions, the Cooper channel 
interactions are always irrelevant. We can see this by noting that such interactions bosonize to 

Vm.m'fjirnflmfLm'fRm' ~ ^ . (115) 

When Of,tot is pinned, the fluctuations of the </> fields are increased compared to the free fermion case and the scaling dimensions 
of the Vm,m' terms are increased. This is a quasi-ID manifestation of the suppression of the Cooper channel by the Amperean 
interaction effects in the gauge theory: oppositely oriented currents repel, hence two oppositely moving fermions in a Cooper 
pair repel and do not like to form such a pair. This stability extends also over a finite range of forward scattering interactions, as 
long as all Vm,m' terms remain irrelevant. 
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The second observation is that if we start with a particle-hole-symmetric and stable fixed point, it means that the scaling 
dimensions of operators are all greater than 2. The fact that in the particle-hole-symmetric case Kn,m' have 

definite phases does not affect such stability considerations. Hence we predict that introducing explicit particle-hole symmetry 
breaking will only add irrelevant perturbations that do not destabilize the phase. This is consistent with our numerical DMRG 
observations on the cylinders where we broke the particle-hole symmetry by allowing tunneling to another layer and observed 
that the central charge of the gapless phase remained unchanged. We also expect such stability to hold in the 2D Dirac CFL 
as well. The absence of the particle-hole symmetry microscopically does manifest itself in observables even though the system 
flows to a fixed point with emergent particle-hole symmetry, since there are no microscopic restrictions for observables to pick 
up components of either formally particle-hole-even or odd combinations in the fixed point theory. Nevertheless, the emergence 
of an effective particle-hole symmetry even when one is not present microscopically makes our numerical studies of strictly 
particle-hole-symmetric models only more important, as a more clear way to access the fixed point theory of the general CFL 
phase. 


Appendix J: Exactly solvable model corresponding to two-wire Dirac CFL 

It is instructive to consider a descendant phase with only two wires cutting through the Fermi sea, which can happen for small 
Ly. We do not find this phase in the model with the screened Coulomb interactions used in the main text; instead, the system 
appears to develop a charge density wave for Ly < 8. The two-wire CFL phase does occur for different more short-range 
interactions and in fact was discovered by Bergholtz and Karlhede in Ref. 56. However, their interpretation does not connect 
with the CFL picture, and we also point out the particle-hole symmetry aspect. 

Bergholtz and Karlhede studied the following electronic model, formulated in the same orbital basis as in Appendix B: 

H = '^[VionjUj+i + V2onjnj+2 - V'2i(c]cj+iCj+2c]+3 -f H.c.)] , (Jl) 

3 

with real and positive Vio, V 20 , and V 21 . The above form taken from Ref. 56 uses different conventions from our Eq. (B9), but 
this is not important. Of all the electron interaction terms, this model keeps only nearest- and next-nearest-neighbor repulsions 
plus the simplest four-fermion term that does not reduce to density-density interaction, and is a natural model near so-called thin 
torus limit. Focusing on the half-filled case, Bergholtz and Karlhede observed that if one considers a subspace spanned by all 
configurations where each pair of sites (2 J, 2 J-f 1) has precisely one electron, the above Hamiltonian acts within this subspace. 
They also showed that, in the regime of interest here, the global ground state resides in this subspace. (More precisely, there are 
two exactly degenerate ground states, and the second one is obtained from the first one by a translation by one lattice spacing, see 
also discussion in Apps. B-C.) Working in this subspace, we can associate a spin-1/2 degree of freedom with each pair of sites 
(2 J, 2 J + 1), where 10 and 01 configurations correspond to spin up and down respectively. We can then identify spin operators 
as Sj = n 2 j — 1/2 = 1/2 — n 2 j+i, Sj = cljC 2 j+i, and the Hamiltonian acting in this subspace is simply an XXZ chain, 

H = ^[( 2 L 2 o - V,o)S}S}+, + 1/21(5+57+1 + H.c.)] . (J 2 ) 

J 

Note that J here refers to the pair (2J, 2J + 1) of the original electronic orbitals, so the spacing between the “sites” in the 
spin chain is twice that in the electronic chain. For 2 V 21 > I 2 V 20 ~ ^ 10 1^ the spin chain is in a critical phase with one gapless 
mode. The dominant power law correlations are staggered along the spin chain, (5+57) ~ /\J ~ and 

(SjSj,) ^ (—1)/“/'/| J — J'P®, where g is Luttinger parameter varying between 1 for the XX chain and 1/2 for the XXX 
chain. There is also a non-oscillating ^ 1/\J — J'p contribution to the (5j5},) correlation, whose power is fixed by the 
conservation of total 5^. 

We now translate these results to the electrons. First, we point out that the 5* correlations correspond to the electron den¬ 
sity correlations at qy = 0, while the 5+ correlations correspond to the density correlations at qy — 2tx j Ly, cf. Apps. B-C. 
Furthermore, since each spin site corresponds to two electronic sites, the staggered correlations in the spin chain correspond 
to correlations at wavevector 7r/2 in the electronic chain, in units where the spacing between orbitals in Eq. (Jl) is 1. In units 
used in the main text, this wavevector corresponds to Q = LyjL This immediately allows us to make interpretation in terms 
of a two-wire CEL with wires at ky = 'n j Ly and —TilLy. Indeed, by the Luttinger theorem discussed in the main text, the 
length of each wire must be Q 1/2 = which is precisely the above Q. Eurthermore, our analysis in the main text predicts 

singularities in the density structure factor at [q^, qy) = {Q 112 , 0) (right to left transfer within a wire) and {Q 1 / 2 , 27r/Ly) (right 
to left transfer between the wires), in agreement with the spin chain solution. Also, from the microscopic expression (in chain 
units) (rTj ~ 1/2)6“*'^“=-^ = (1 — we can readily deduce that the electron structure factor at gj, = 0 and 

small qx behaves as in agreement with the prediction at the end of App. 13 for the general quasi-lD case. Einally, the exact 
ground state resides in the sector with the 0110 root configuration,^® which is what we found for all cases with even number of 
wires (APBC). Thus, we can match properties of the exact solution with those expected from the two-wire Dirac-CEL. 
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It is instructive to identify the microscopic symmetries of the electronic problem discussed in App. B with symmetries of the 
above spin chain. This will also lead us to a more microscopic identification of the composite fermions. The anti-unitary mirror 
symmetry Eq. (B8) translates to 

M^T : S} S^j , Sj Sj , SJ SJ , (13) 

and hence acts like a “boson time reversal” when the spin model is interpreted as a hard-core boson model. 

For the spatial inversion, we consider inversion in the bond center (midpoint between two neighboring orbitals), Eq. (B7). 
This preserves the above restricted Hilbert space and becomes in the spin model 

I' : S^j -Slj , Sj SZj , SJ Sjj . (14) 

Most interestingly, the anti-unitary particle-hole symmetry Eq. (A3) becomes 

m : S^j -S} , Sj -SJ , SJ -Sj , (15) 

i.e., T’H : Sj —> —Sj, which acts like familiar spin time reversal. 

We can now guess what the composite fermions are in the quasi-lD descendant Dirac-CFL theory in this case. They can 
be viewed as arising from a slave particle description of the spin chain, where we write Sj = with a constraint 

+ flifjl = 1 on each site, and postulate a mean field where the spinons and hop independently. Since Sj 
carries transverse momentum 27r/Lj,, we can interpret as carrying momentum ky = n/Ly and as carrying -■K/Ly, similar 
to the two-wire picture with APBC in the transverse direction. Furthermore, we can implement all of the above symmetries 
in terms of the spinons and match these with the proposed transformation properties in the quasi-ID Dirac-CFL in App. I. 
In particular, we can implement the anti-unitary particle-hole symmetry in the familiar way as time reversal symmetry on the 
spinons, fj^ fj^, fj^ —/jt> which immediately gives us the Kramers doublet physics. [More precisely, one needs to 
modify the spinon right- and left-moving fields with appropriate phase factors to match with the transformation properties of 
the Dirac-CFL fields. When matching, one needs to remember that the inversion symmetry I' in this Appendix is relative to a 
mid-point between two orbitals, while in App. I the inversion I is relative to the origin which coincides with one orbital; I and 
/' are related by translation by one orbital.] 

Interestingly, since microscopically Sj = c\jC 2 j+i and we are considering subspace with cljC 2 j + c\j_^^C 2 j+i = 1, we 
can interpret the spinons as the original fermions in the ELL orbitals on the cylinder that develop spontaneous coherence along 
the chain independently on the even and odd sublattices, writing schematically fj^ ^ C 2 j, fji ~ C 2 j+i. Note that while 
microscopically the T’H maps C 2 j —?> c\ j and C 2 j+i —in the considered subspace this action coincides with the above 
familiar time reversal transformation of the spinons. While the above is special to the two-wire case, it would be interesting to 
find a microscopic derivation of the quasi-ID Dirac-CFL fields also for cases with more wires. 

We can also examine what happens when we break the particle-hole symmetry and the interplay with the other symmetries. 
For example, we can consider a perturbation 


6H = veiclcj+lCj+ 2 Cj _^_3 + H.c.)(l - nj^i - nj+ 4 ) , (J6) 

3 

which preserves the anti-unitary mirror and the inversion symmetries but breaks the particle-hole symmetry. This perturbation 
also acts in the above restricted Hilbert space and becomes in the spin model 

SH = J2veiSjSj^, + H.c.)( 5}+2 - ^5-i) ■ (J7) 

3 

We can understand its effects using bosonization. Let us use a hydrodynamic description of the XXZ chain where the S^ spin 
component is represented as Sj = dxO + A{—lY sin(20) and the bond energy [J, J -|- 1] is represented as B{—\Y cos(20). 
In these variables, additional S^S’^ interactions in the XXZ chain that can gap it out (a.k.a. “umklapp” terms in half-filled boson 
or fermion chains) contribute Acos(40); this term of course respects all the symmetries of the XXZ chain, but in the critical 
phase it must be irrelevant. Now we see that when we break the particle-hole symmetry by the above 5H, it contributes a term 
A' sin(46*), which must also be irrelevant as long as the unperturbed critical phase is stable. This observation connects with our 
discussion of particle-hole symmetry conditions on for four-fermion terms in the general Dirac-CFL theory in App. 14 (with 
different bosonization conventions between the two sections, since in App. 14 we effectively grouped right spin-up with left 
spin-down fields). 

As another example, we can consider a perturbation 

5H' = ^r;^(c]cj+iCj+ 2 cJ +3 -f H.c.)(l/2-nj_i) = - ^ -f H.c.)S'}_i . (18) 

j j 


27 


This perturbation preserves the anti-unitary mirror symmetry but breaks both the inversion and particle-hole symmetry. In the 
spin model, we see that the effect of this perturbation is roughly similar to adding a field in the direction on the XXZ spin 
chain. The effect on the correlations is that the 5'+ correlations remain staggered, while the correlations are shifted to an 
incommensurate wavevector. This is consistent with developing an imbalance between the spin-up and spin-down spinons— 
correspondingly, between the ky = 'k jLy and —tt/Lj, wires, which is expected when both the inversion and particle-hole 
symmetries are broken. We thus see that the Bergholtz and Karlhede model realizing the simplest quasi-ID descendant of the 
CFL provides useful playground for examining the interplay of symmetries in the half-filled Landau level problem. 



